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We argue that screening of higher-representation color charges by gluons implies a domain structure in the 
vacuum state of non-abelian gauge theories, with the color magnetic flux in each domain quantized in units 
corresponding to the gauge group center. Casimir scaling of string tensions at intermediate distances results 
from random spatial variations in the color magnetic flux within each domain. The exceptional G(2) gauge 
group is an example rather than an exception to this picture, although for G(2) there is only one type of vacuum 
domain, corresponding to the single element of the gauge group center. We present some numerical results for 
G(2) intermediate string tensions and Polyakov lines, as well as results for certain gauge-dependent projected 
quantities. In this context, we discuss critically the idea of projecting link variables to a subgroup of the gauge 
group. It is argued that such projections are useful only when the representation-dependence of the string 
tension, at some distance scale, is given by the representation of the subgroup. 



I. INTRODUCTION 

Over the past decade, a great deal of numerical evidence has 
been obtained in favor of the center vortex confinement mech- 
anism (for reviews, cf. refs. [1] and [2]). According to this 
proposal, the asymptotic string tension of a pure non-abelian 
gauge theory results from random fluctuations in the number 
of center vortices linked to Wilson loops. It is sometimes 
claimed that gauge theory based on the exceptional group 
G(2) is a counterexample to the vortex mechanism, in that 
the theory based on G(2) demonstrates the possibility of hav- 
ing "confinement without the center" [3]. This claim has two 
questionable elements. In the first place, the asymptotic string 
tension of G(2) gauge theory is zero, in perfect accord with 
the vortex proposal. No center vortices means no asymptotic 
string tension. From this point of view G(2) gauge theory is 
an example of, rather than a counterexample to, the impor- 
tance of vortices. The question of whether G(2) gauge theory 
is nonetheless confining then depends on what one means by 
the word "confinement". This is a basically a semantic issue, 
but for the sake of clarity we will explain our view below. The 
second point is that G(2), like any group, does have a center 
subgroup. This may seem like a quibble, since the center of 
G(2) is trivial, consisting only of the identity element. We 
will argue, however, that this group element has an important 
role to play, not only in G(2) but perhaps also in SU(N) gauge 
theories, in the classification of non-trival vacuum structures. 

We begin with a discussion of the term "confinement", 
which is used in the literature in several inequivalent ways: 

1 . Confinement refers to electric flux-tube formation, and 
a linear static quark potential. 

2. Confinement means the absence of color-electrically 
charged particle states in the spectrum. 

3. Confinement is the existence of a mass gap in a gauge 
field theory. 



We prefer the first of these options. Real QCD, according to 
this definition, confines at intermediate but not at large dis- 
tance scales, since electric flux tubes break, the static poten- 
tial goes flat, and excited hadronic states with string-like con- 
figurations are only metastable. For this reason we use the 
term "temporary confinement" to describe the situation in real 
QCD with light quarks. Most discussions, however, do not 
make such a distinction, and apply the label "confinement" 
to real QCD without any qualification. In that case, to con- 
form to common usage, why not adopt the second or third 
definitions? Our reason is that this choice would force us to 
also describe many other theories as "confining" which are not 
normally considered so, and which (unlike real QCD) have 
nothing at all to do with flux tube formation and confining 
potentials. 

As one example, let us consider an electric superconductor, 
or, in its relativistic version, the abelian Higgs model. In this 
theory, any external electric charge is screened immediately 
by the condensate. There is a mass gap, and asymptotic par- 
ticle states are all electrically neutral. If definitions (2) and/or 
(3) define confinement, then the abelian Higgs model is surely 
confining. But this contradicts standard usage in the literature, 
where confinement of electric charge is often identified with a 
dual (i.e. magnetic) abelian Higgs model, while the ordinary 
abelian Higgs model is identified with confinement of mag- 
netic, not electric, charge. 

A second example along the same lines is an electrically 
charged plasma. In response to an external charge, the plasma 
rearranges itself to screen out the field (and hence the charge) 
of the probe. Since all probe charges are neutralized, and 
all Coulombic fields are Debye screened, an electrical plasma 
could also be taken, according to definition (2), as a confining 
system. 

Our third example is an SU(2) gauge-Higgs theory with a 
pair of Higgs fields in the adjoint representation. In this case 
there are two distinct massive phases of the system, charac- 
terized by an asymptotic string tension which is non-zero in 
one phase, and zero in the other. We have every right, in this 
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theory, to distinguish between the Confinement phase and the 
Higgs phase. But definitions (2) and (3) allow no such dis- 
tinction. Both phases, according to these definitions (and ac- 
cording to the Fredenhagen-Marcu criterion [4]) are equally 
confining. 1 Again, we feel that this is an abuse of the term 
"confinement". 

Our point is that consistent application of any of the three 
suggested definitions does at least some violence to common 
usage. According to our preferred terminology, i.e. option (1), 
both real QCD and G(2) gauge theory are examples of only 
"temporary" confinement. The other options, in our opinion, 
are worse. We do not think that electric superconductors, elec- 
trical plasmas, and gauge systems which are clearly in a Higgs 
phase, should be described by a term which is so often and so 
generally associated with electric flux tubes and a linear po- 
tential. 

Of course, the serious question to be addressed is not 
whether or not to refer to G(2) gauge theory as confining; at 
the end of the day this is only a matter of words. The real 
challenge which is posed by G(2) gauge theory, in our opin- 
ion, is the fact that this theory has a non-vanishing string ten- 
sion for the static quark potential over some finite range of 
distances, intermediate between the breakdown of perturba- 
tion theory and the onset of color screening. Where does that 
string tension come from, if not from vortices? This ques- 
tion has, in fact, a much broader context. In connection with 
SU(N) gauge theories, we have emphasized many times in the 
past (beginning with ref. [5]) that the dependence of the string 
tension on the color group representation of the static quarks 
falls into two general categories, depending on distance scale: 

1. The Casimir Scaling Region — At intermediate dis- 
tances, the string tension of the static quark potential, 
for quarks in group representation r, is approximately 
proportional to the quadratic Casimir C r of the repre- 
sentation, i.e. 

C r 

O r ~—0 F (1.1) 
Lf 

where F denotes the fundamental, defining represen- 
tation. Thus, e.g., the string tension of quarks in the 
adjoint representation of SU(N) gauge theory is about 
twice as large as that of quarks in the fundamental rep- 
resentation. 

2. The Asymptotic Region — Asymptotically, due to color 
screening by gluons, the string tension depends only on 
the transformation properties of representation r with 
respect to the center subgroup; i.e. on the N-ality k r in 
an SU(N) gauge theory: 

o r ^a(k r ) (1.2) 



where a(k) is the string tension of the lowest dimen- 
sional representation of SU(N) with N-ality k. These 
are known as the "k-string tensions." The string tension 
of adjoint representation quarks in the asymptotic re- 
gion is zero. 2 

We have no reason to think that G(2) gauge theory dif- 
fers from SU(N) gauge theories with respect to the interme- 
diate/asymptotic behaviors of the string tension. First of all, 
there is certainly a linear static quark potential in an interme- 
diate region (as we confirm by numerical simulations), and 
we think it likely that the representation dependence of the 
string tension at intermediate distances follows an approx- 
imate Casimir scaling law. Casimir scaling has not been 
checked numerically for G(2), but for the purposes of this 
article we will assume it to be true, pending future investi- 
gation. Secondly, all G(2) group representations transform in 
the same way with respect to the center subgroup, i.e. triv- 
ially, and all must have the same asymptotic string tension, 
namely zero, because of color screening by gluons. There- 
fore, the existence of a linear potential in G(2) is the subject 
of a broader question, relevant to any gauge group: If center 
vortices (or the absence of center vortices) explains the values 
of the asymptotic string tension, then what accounts for the 
linear static quark potential at intermediate distance scales, 
where the string tension depends on the quadratic Casimir of 
the gauge group, rather than the representation of the gauge 
group center? 

In the next two sections we will present our answer to this 
question, which improves, in some crucial ways, on an old 
model introduced by two of the present authors, in collab- 
oration with M. Faber [6]. We begin by asking (section II) 
what features, of a typical gauge field vacuum fluctuation, can 
possibly account for both Casimir scaling and color screen- 
ing. We are led to postulate (in section III) a kind of domain 
structure in the vacuum, with magnetic flux in each domain 
quantized in units corresponding to the elements of the cen- 
ter subgroup, including the identity element. In section IV we 
present some lattice Monte Carlo results for G(2) gauge the- 
ory, which confirm (as expected) the existence of a string ten- 
sion over some intermediate distance range. In section V we 
consider fixing the gauge in G(2) lattice gauge theory, so as to 
leave either an SU(3), Z3, or SU(2) subgroup unfixed, together 
with a projection of the G(2) link variables to the associated 
subgroup. While there are some numerical successes with this 
procedure, suggestive of SU(3) or Z3 or SU(2) "dominance" 
in G(2) gauge theory, in the end we were unable, at least thus 
far, to find any correlation between the projected observables 
and gauge-invariant observables. We think this failure makes 
a point which may also be relevant to, e.g., abelian projec- 



The Fredenhagen-Marcu criterion, also advocated by the authors in ref. 
[3] as a confinement criterion, really only tests whether the charge of an 
external source is screened by dynamical matter fields, an effect which 
also occurs in a plasma and an electric superconductor. 



2 Of course, it may be that the k-string tensions are proportional to C rmin , 
where r mm is the lowest dimensional representation of N-ality k. Con- 
fusingly, this hypothesis is also called "Casimir scaling" in the literature. 
However, the term "Casimir scaling" was originally introduced in ref. [5] 
to refer to string tensions in the intermediate region, obeying eq. (1.1), and 
we will stick to that original definition here. 
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tion in SU(N) theories: Projection to a subgroup is unlikely 
to be of much significance, unless the group representation- 
dependence of the string tension, at some distance scale, de- 
pends only on the representation of the subgroup. Section VI 
contains some concluding remarks. 



II. CASIMIR SCALING 

Casimir scaling of string tensions is quite natural in non- 
abelian gauge theory. Consider a planar Wilson loop in group 
representation r. Take the loop to lie in, e.g., the x — y plane at 
z = t = 0, and imagine integrating over all links, in the func- 
tional integral, which do not lie in the plane, i.e. 

W r (C) = ^ J DU,(x,y,z,t)x r [U(C)]e- s * 

' ' DU x (x,y,0,0)DU y (x,y,0,0) % r [U(C)} 

[DU x DU y ] z and/or t ^ Q 
DU z {x,y,z,t)DU t (x,y,z,t) e~ Sw } 
= l - J DU x (x,y)DU y (x,y) % r [U{C)] 



x exp 



S e ff[U x ,Uy, 



(2.1) 



where % r is the group character in representation r, and U (C) 
the loop holonomy. Then the computation of the planar loop 
in D = 4 dimensions, with the Wilson action Sw, is reduced to 
a computation in D = 2 dimensions, with an effective action 
S e ff. This effective action can be regarded as the lattice regu- 
lation of some gauge-invariant D = 2 dimensional continuum 
action S c f^. Although S c °^ is surely non-local, non-locality 
can, in some circumstances, be traded for a derivative expan- 
sion, so we have 

S c e ° f n f = f d 2 x (a Q Tr[F 2 }+a l Tr[D p FD fl F}+a 2 Tr[(F 4 } + ...^ 

(2.2) 

where F = F xy . Truncating S L e 'jj to the first term, which ought 
to be dominant at large scales, we obtain 

W r (C) ~ X - J DA x (x,y)DA y (x,y) % r [U(C)} 



x exp 



-ClQ 



I d 2 x Tr[F 2 



(2.3) 



Since this is simply gauge theory in D = 2 dimensions, the 
functional integral can be evaluated analytically, and the re- 
sult is that the string tension is proportional to the quadratic 
Casimir C r . It should be kept in mind, however, that the va- 
lidity of the derivative expansion in (2.2) is an assumption. 
There might be some non-local terms in S e ff which are not 
well represented by such an expansion. 

The argument for an effective D — 2 action can also be 
phrased in terms of the ground state wavefunctional ^^[A] 
of D+ 1 dimensional gauge theory in temporal gauge. It was 



originally suggested in ref. [7] that for long-wavelength fluc- 
tations, the ground state could be approximated as 



<[A] 



exp 



/ j d D x Tr[F 2 



(2.4) 



In that case 
W r (C) 



{ Xr [U(C)}) D =^^l\Xr[U(C)}\^) 

(X r {U(C)}) D = 3 = (V 2 \x r {U(C)]\V 2 ) 

(Xr[U(C)}} 



D=2 



(2.5) 



and again the D — 4 dimensional calculation has been reduced, 
in two steps, to the D = 2 dimensional theory, where the string 
tensions are known to obey the Casimir scaling law. 3 The ar- 
gument hangs on the validity of the approximation (2.4). This 
approximation is supported by Hamiltonian strong-coupling 
lattice gauge theory, where the ground state can be calcu- 
lated analytically [13], and it has also been checked numer- 
ically [14]. Karabali, Kim, and Nair [15] have pioneered an 
approach in which the ground state wavefunctional ^[A] of 
D = 2 + l dimensional gauge theory can be calculated analyt- 
ically in the continuum, as on the lattice, in powers of 1 /g 2 , 
and according to these authors the leading term is again given 
by eq. (2.4). See also Leigh et al. [16] and Freidel [17]. 

Returning to eq. (2.3), we note that what it really says is 
that on a D = 2 slice of D = 4 dimensions, the field strength 
of vacuum fluctuations is uncorrected from one point to the 
next. This is because, in D = 2 dimensions, there is no Bianchi 
identity constraining the field strengths, and, by going to an 
axial gauge, the integration over gauge potentials A x , A y can 
be traded for an integration over the field strength F^. The 
absence of derivatives of F in eq. (2.3) means that the field 
strengths fluctuate independently from point to point. This 
is, of course, absurd at short distance scales, and results from 
dropping the derivative terms in eq. (2.2). On the other hand, 
suppose that Fxy has only finite-range correlations, with cor- 
relation length I on the two-dimensional surface slice. In that 
case, dropping the derivative terms may be a reasonable ap- 
proximation to S e ff at scales larger than /, and Casimir scaling 
is obtained at large distances. 

Regardless of the derivation, the essential point is the fol- 
lowing: Casimir scaling is obtained if field strengths on a 2D 
slice of the 4D volume have only finite range correlations. 4 
Then the color magnetic fluxes in neighboring regions of area 
I 2 are uncorrected, the effective long-range theory is D=2 di- 
mensional gauge theory, and Casimir scaling is the natural 
consequence. 

Once Casimir scaling gets started, the next question is why 
it ever ends; i.e. why the representation dependence suddenly 
switches from Casimir scaling to N-ality dependence. The 



3 Dimensional reduction from four to two dimensions was advocated inde- 
pendently, on quite different grounds, by Olesen [8], who (together with 
Ambj0rn and Peterson [9]), noted that this reduction implies eq. (1.1). See 
also Halpern [10] and Mansfield [11]. 

4 In this connection, see also Shoshi et al. [12]. 
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conventional answer is based on energetics: As the flux tube 
gets longer, at some stage it is energetically favorable to pair 
create gluons from the vacuum. These gluons bind to each 
quark, and screen the quark color charge from representation 
r to representation r mi „ with the lowest dimensionality and 
the same N-ality as r. In particular, if r m! „ is a singlet, then 
the flux tube breaks. Note that this means that the approxi- 
mations (2.3) and (2.4), which imply Casimir scaling, must 
break down at sufficiently large distance scales (at least for 

Ncolors < °°)- 5 

While the energetics reasoning for color screening via glu- 
ons is surely correct, when we speak of the "pair creation" of 
gluons out the vacuum, we are using the language of Feyn- 
man diagrams, and describing, in terms of particle excita- 
tions, the response of the vacuum to a charged source. How- 
ever, the path-integral itself is a sum over field configurations. 
When a Wilson loop at a given location is evaluated by Monte 
Carlo methods, the lattice configurations which are generated 
stochastically have no knowledge of the location of the Wil- 
son loop being measured, and no clue of where gluon pair 
creation should take place. The Wilson loop is just treated 
as an observable probing typical vacuum fluctuations, not as 
a source term in the action. Nevertheless, the Monte Carlo 
evaluation must somehow arrive at the same answer that is 
deduced, on the basis of energetics, from the particle picture 
of string breaking. This leads us to ask the following ques- 
tion: What feature of a typical vacuum fluctuation can account 
for both Casimir scaling of string tensions at intermediate dis- 
tances, and center-dependence at large distances? What we 
are asking for is a "field" explanation of screening, to com- 
plement the particle picture. 

An important observation is that if the asymptotic string 
tension is to depend only on N-ality, then the response of a 
large Wilson loop to confining vacuum configurations must 
also depend only on the N-ality of the loop. Center vortices 
are the only known example of vacuum configurations which 
have this property. The confinement scenario associated with 
center vortices is very well known: random fluctuations in 
the number of vortices topologically linked to a Wilson loop 
explain both the existence and the N-ality dependence of the 
asymptotic string tension. Casimir scaling, on the other hand, 
seems to require short-range field strength correlations with 
no particular reference to the center subgroup. As Casimir 
scaling and N-ality dependence each have straightforward, but 
different explanations, the problem is to understand how both 
of these properties can emerge, in a lattice Monte Carlo simu- 
lation, from the same set of thermalized lattices. 



HI. A MODEL OF VACUUM STRUCTURE 

Since Casimir scaling is an intermediate range phe- 
nomenon, while N-ality dependence is asymptotic, we suggest 



5 For a treatment of the long-range Yang-Mills effective action in the context 
of strong-coupling lattice gauge theory, cf. ref. [18]. 



that the random spatial variations in field strength, required for 
Casimir scaling, are found in the interior of center vortices. 
These interior variations cannot be entirely random however; 
they are subject to the constraint that the total magnetic flux, 
as measured by loop holonomies, corresponds to an element 
of the gauge group center. 

The picture we have in mind is indicated schematically in 
Figs. 1 and 2; these are cartoons of the vacuum in a 2D slice 
of the D = 4 spacetime volume. The first figure is an im- 
pression of the vacuum in Yang-Mills theory, the second is 
for G(2) theory. Each circular domain is associated with a to- 
tal flux corresponding to an element of the center subgroup. 
The domains may overlap, with the fields in the overlap re- 
gions being a superposition of the fields associated with each 
overlapping domain. In the SU(2) case there are two types of 
regions, corresponding to the two elements of the center sub- 
group. Regions associated with the Z2 center element z = — 1 
are 2D slices of the usual thick center vortices, and the straight 
lines are projections to the 2D surface of the Dirac 3 -volume 
associated with each vortex. But domains corresponding to 
the z — 1 center element are also allowed; we see no reason 
to forbid them. Within each domain, the field strength fluc- 
tuates almost independently in sub-regions of area I 2 , where 
/ is a correlation length for field-strength fluctuations, subject 
only to the weak constraint that the total magnetic flux in each 
domain corresponds to a center element of the gauge group. 
In the case of G(2) there is of course only one center element 
Z = 1 of the trivial Z\ subgroup, and no Dirac volume, hence 
the difference in the two cartoons. This picture is sufficient to 
explain the N-ality dependence of loops which are large com- 
pared to the domain size, and Casimir scaling of loops which 
are small compared to the domain size. 




FIG. 1 : A 2D slice of the D=4 Yang-Mills vacuum. Circular regions 
with (Dirac) lines correspond to z = — 1 domains, circular regions 
without lines denote z = +1 domains. 

To support the last statement we return to an old model in- 
troduced in ref. [6]. In this model, it is assumed that the effect 
of a domain (2D cross-section of a vortex) on a planar Wil- 
son loop holonomy is to multiply the holonomy by a group 
element 

G(a n ,S)=Sexp[ia n -H]S f (3.1) 

where the {Hi} are generators of the Cartan subalgebra, S is 
a random group element, oc" depends on the location of the 
domain relative to the loop, and n indicates the domain type. 
If the domain lies entirely within the planar area enclosed by 
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FIG. 2: A 2D slice of the D=4 vacuum of G(2) gauge theory. There 
is only one type of domain, corresponding to the single element of 
the center subgroup. 



the loop, then 



where 



exp[id n -H) =zj 



(3.2) 



If we make this assumption, then 



(Wj(C)) 



n^ 1 +fog J [<%(x)]}wf ert (C) 



= exp 



£tog{(l-/i-/o) + /i^[o£(*)] 

. X 

+f Q gj[a? c (x)]}]wf ert (C) 



(3.6) 



where the product and sum over x runs over all plaquette posi- 
tions in the plane of the loop C, and a£(x) depends on the po- 
sition of the vortex midpoint x relative to the location of loop 
C. The expression Wj en '(C) contains the short-distance, per- 
turbative contribution to Wj (C) ; this will just have a perimeter- 
law falloff. 

To get the static potential, we consider loop C to be a rect- 
angular RxT contour with T S> R. Then the contribution of 
the domains to the static potential is 



Zn = e 2*in/N eZN 



(3.3) v j (R) = - 1 £ log{(l-/i-/o)+/i^[ai(x)]+/o^[ag(x)]} 



and / is the unit matrix. At the other extreme, if the domain is 
entirely outside the planar area enclosed by the loop, then 



exp[id"-H]=I 



(3.4) 



For a Wilson loop in representation r, the average contribution 
from a domain will be 



g r [a n ]Id r = JdSSexp\ia n -H]S f 



d r 



exp [/a" H] 



(3.5) 



where d r is the dimension of representation r, and Id r is the 
unit matrix. 

Consider, e.g., SU(2) lattice gauge theory, choosing H = 
L3. The center subgroup is Z2, and there are two types of 
domains, corresponding to zn = 1 and z\ = — 1. Let f\ repre- 
sent the probability that the midpoint of a z\ domain is located 
at any given plaquette in the plane of the loop, with fo the 
corresponding probability for a zo = 1 domain. 6 Let us also 
make the drastic assumption that the probabilities of finding 
domains of either type centered at any two plaquettes x and y 
are independent. Obviously this "non-interaction" of domains 
is a huge over-simplification, but refinements can come later. 



6 In the continuum, /o, 1 become probability densities to find the midpoints 
of domains at any given location. These should not be confused with the 
fraction of the lattice covered by domains, which depends also on domain 
size. It is possible that every site on the lattice belongs to some domain, or 
to more than one overlapping domains. 



(3.7) 

where x is now the distance from the middle of the vortex to 
one of the time-like sides of the Wilson loop. The question is 
what to use as an ansatz for a^(x). One possibility was sug- 
gested in ref. [6], but there the desired properties of the poten- 
tial in the intermediate region, namely linearity and Casimir 
scaling, were approximate at best. We would like to improve 
on that old proposal. 

Our improved ansatz is motivated by the idea that the mag- 
netic flux in the interior of vortex domains fluctuates almost 
independently, in subregions of extension /, apart from the 
restriction that the total flux results in a center element. To 
estimate the effect of this restriction, consider a set of random 
numbers {F„},n = 1,2, ...,M, whose probability distributions 
are independent apart from the condition that their sum must 
equal 2tc, i.e. 



«2[{Fn}}) 



lUndF n Q[{F n }}e-^ F ^ (Lf =1 F n - 2k 



mndF n e-^i Hl M^ Fn _ 2% 



Then for a sum of A < M random numbers 

A A 2 
M~M 2 




271 — 

M 



(3.8) 



(3.9) 



If the constraint is, instead, that £ n=1 F„ = 0, then the sec- 
ond term on the rhs is dropped. We then postulate that the 
OCc(x) phase likewise consists of a sum of contributions from 
subregions in the domain of area I 2 , which lie in the interior 
of loop C. It is assumed that these contributions are quasi- 
independent in the same sense as the {F n }. Then if the total 
cross-sectional area of a vortex overlapping the minimal area 
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of loop C is A 3> I , we take as our ansatz 



(«£(*)) 



2a, 



A A 2 



271- 



A v 



(3.10) 



where A V ,A' V are the cross-sectional areas of the n = 1 and 
n = domains, respectively. For the sake of simplicity, let 
us suppose that A' v = A v , with / = 1 lattice spacing, and also 
imagine that the cross-section of a vortex is an L v x L v square. 
In that case 



Vj(R) = £ log{(l-/i-/o)+/i^[aiWl 
l 



+/o£y [<#(*)]} 



(3.11) 



Now consider two limits: small loops with R <C L v , and large 
loops with R » L v . In the former case 



2^v 

Vj(R) « - £ log{(l-/i-/o) + /i^-[ai(x)] 

+/o£y [<$(*)]} 
« -ivlog{(l-/i-/ )+/i^[ai] 

+fogjK]} 



(3.12) 



where 



(4) : 



But for small 9 



Z 2 f R R 2 
2^ 



e 2 



2tc— 



(3.13) 



^■(e)»i-g-7C/ + i) 



(3.14) 

so putting it all together, we get for R^L v a linear potential 

Vj(R) = ^A^j(j+W (3.15) 
with a string tension 

f, 4- fr. T 2 

(3.16) 



6 2,1, 



which is proportional to the SU(2) Casimir. In the opposite 
limit, for R > L v , 



Vj(R) 



£ log{ (l - /i - /o) + /i gj[a l R ] + fyg } [<\ } 



-- 2 Lv 



-J?log{(l-/i)+/i^[27i]} 



(3.17) 



where the summation runs over domains which lie entirely 
within the minimal loop area, so that = 2%, = 0. For 
those values, £/(oc^) = — 1 for j = half-integer, £? 7 (a^) = 1 
for j = integer, and £/(oc^) = 1 for all j. Then Vj is again a 
linear potential, with asymptotic string tension 



asy _ 



■ log ( 1 - 2/i ) j = half-integer 
j = integer 



(3.18) 



This has the correct N-ality dependence. One can adjust the 
free parameter /j to get a potential for the j =\ representation 
which is approximately linear, with the same slope, at all R. 

In generalizing to SU(N) we must take into account that 
there are N types of center domains, so that 



(Wr(C)} 



N-l 



EI l-I/n(l-Re^W]) \wr rt (c) 

x I n=0 I 



= exp 



where 



iV-l 



DogU-L/„(l-Re^(x)]) 

x \ n=0 



exp [/a" H] 



and 



fn=f N -n , QrW} = Q;W- n } 



W r pert (C) 

(3.19) 
(3.20) 
(3.21) 



For very small loops, which are associated with small a, we 
have that 



1 



l-g r [a n ] w ±a?a n j-Tr[H'HJ] 

U r 

1 



2(JV 2 - 1) 



a" • a"C r 



(3.22) 



where C r is the quadratic Casimir for representation r. Once 
again, if d" • d" is proportional to the area of the vortex in the 
interior of the loop, for small loops, we end up with a linear 
potential whose string tension is proportional to the quadratic 
Casimir. For large loops, if the vortex domain lies entirely in 
the loop interior, we must have 



1 



exp[/d" H] 



(3.23) 
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where k is the N-ality of representation r. Again we obtain a 
linear potential, whose string tension depends only on the N- 
ality. It is possible to choose the /„ so that for very large loops 
the k-string tensions follow either the Sine Law, or are propor- 
tional to the Casimir C Tmin ; c.f. ref. [19] for a full discussion of 
this point. 

Essentially, we are proposing a model of vacuum structure 
which is much like an instanton or monopole or caloron gas 
picture. The essence of such models is that the functional in- 
tegral over all configurations is approximated, in the infrared, 
by summation over a special set of configurations, and fluc- 
tuations around those configurations. In our case, the special 
set is center domains. A center domain, whether in SU(N) or 
G(2) gauge theory, is a planar region and a constraint. The 
constraint is that color magnetic fluctuations within the pla- 
nar region add up to a center element, as defined by a loop 
holonomy around the boundary of the domain. The vacuum 
fluctuations of a gauge field in a plane is then approximated 
by a sum over domain positions, and an integral over the (con- 
strained) magnetic fluctuations within each domain. This idea 
can surely be elaborated and perfected; the implementation in 
this section is only a first, and certainly over-simplified, at- 
tempt at its realization. 

To illustrate how things may go in SU(2), we choose the 
following parameters 

• L v = 100 (square cross section) 

• f x = 0.01, f = 0.03 
. ^/(2/i)=4 

These parameters are chosen as an example; we are not sug- 
gesting that they are necessarily realistic. We are also taking 
/ = 1, so the potentials are linear from the beginning. Fig. 3 
shows the resulting potentials for Vj(R) with j = ^,1, | in the 
range R e [0,200]; one clearly sees the N-ality dependence at 
large R. Fig. 4 shows the same data in a region of smaller 
R < 50; there does appear to be an interval whether the po- 
tential is both linear and Casimir scaling. For fixed /o,/i the 
fj parameter is fine-tuned so that the string tension at small 
and large R are about the same, and the potential for j = ^ is 
roughly linear over the entire range of/?. 7 

Finally we turn to G(2). Once again, the logic is that 
Casimir scaling + screening imply a domain structure. The 
model is much like SU(2) except there is only one type of 
domain, corresponding to the single element of the Z\ center 
subgroup. The requirement that the asymptotic string tension 
vanishes for all representations requires that a = when the 
domain lies entirely in the loop interior, so that 



f=.01,g=.03,p=2 



ococ= — 
2p 



A 



Al 



(3.24) 



7 The string tension shown in Figs. 3 and 4 is linear even at the smallest 
values of R because we have taken the field-strength correlation length to 
be / = 1 lattice spacing, and dropped the W pert (C) contribution, which 
contains the perturbative contributions to the static potential. 




200 



FIG. 3: The static potential for j = j, 1, § static sources, for vortex 
width = 100, in the distance range R € [0,200]. 



f=.01,g=.03,p=2 



> 




FIG. 4: Same as Fig. 3, in the range R e [0,50]. 



Together with the fact that Ti r [H l H^] °c 8,jC, in G(2) also, we 
obtain a linear potential with Casimir scaling at small R, and a 
vanishing string tension for large loops, in all representations, 
at large R. 

We note again that Casimir scaling, at the outset of con- 
finement, has never been verified in G(2) lattice gauge theory. 
However, an initial interval of Casimir scaling (followed by 
color screening) seems to be inescapable in our model, and 
can be treated as a prediction. If this scaling is not eventually 
confirmed by lattice simulations in G(2) pure gauge theory, 
then the model is wrong. 

The last question is how the zq = 1 domains on a 2D slice 
extend into the remaining lattice dimensions, in either G(2) 
or SU(N) theories. The z n 7^ 1 domains, in SU(N) theories, 
extend to line-like objects in D = 3 dimension, or surface- 
like objects in D = 4 dimensions, bounding a Dirac surface 
or volume, respectively. These are the usual (thick) center 
vortices. We think it likely that zq = 1 domains also extend 
to line-like or surface-like objects in D = 3 and 4 dimensions. 
This opinion is based on the fact that the Euclidean action of 
vortex solutions, for both z 7^ 1 and z = 1, has been shown to 
be stationary at one loop level [20-22]. While the stability of 
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these solutions, as well as the validity of perturbation theory 
at the relevant distance scales, have not been fully established, 
the existing results do suggest that the z = 1 vortices are not 
too different in structure from their z ^ 1 cousins. 

To summarize: We have improved on the model introduced 
in ref. [6] in two ways. First, we have motivated an ansatz for 
the Oc£.(x) which gives precise linearity and Casimir scaling for 
the static quark potential at intermediate distance scales ("in- 
termediate" begins at R = I). Secondly, we have allowed for 
the existence of domains corresponding to the center element 
Zo = 1 ; this addition permits a unified treatment of SU(N) and 
G(2) gauge theories. The strength of our improved model 
is that it gives a simple account of the linearity of the static 
potential in both the intermediate and asymptotic distance 
regimes, with Casimir scaling in the former and N-ality de- 
pendence in the latter. Its main weakness, in our view, is that 
there is no explanation for why the string tension of funda- 
mental representation sources should be the same at interme- 
diate and asymptotic distances. This equality can be achieved, 
but not explained, by fine-tuning one of the parameters in the 
model. 



IV. STATIC PROPERTIES OF G(2) GAUGE THEORY 

In this section we will present some numerical results con- 
cerning the intermediate string tension and Polyakov line val- 
ues of G(2) lattice gauge theory, determined via Monte Carlo 
simulations. 



At each link we choose A = Z)jfc(oc) for k = 1 . . .7. Define 



A. Metropolis Algorithm 



The partition function is 



= f DUp expJ | £ t 



(4.1) 



where U^(x) are G(2) group elements, and P^(x) is the corre- 
sponding plaquette variable. In a Metropolis Update step, we 
monitor the change of action under a change of a particular 
link: 



Ufa) = A(a) U^x) , 



(4.2) 



where A(a) € G(2) is given in terms of the Euler parameteri- 
zation (B3) as function of the 14 Euler angles, denoted collec- 
tively by d. These angles are chosen stochastically, in a range 
which is tuned to achieve an acceptance rate of roughly 50%. 

In order to reduce auto-correlations as well as to reduce the 
number of Monte-Carlo steps, so-called "micro-canonical re- 
flections" are useful: After a certain number of Metropolis 
steps, the lattice configurations is replaced by a configuration 
of equal action. This is done by a sweep through the lattice, 
and at each link making the replacement 



(4.4) 



where the sum over positive and negative v runs over all 
plaquettes containing link (x,fi). Then the replacement link 
changes the action by an amount 



AS = 



P 



trA(a)B^(x) - tvBfj(x) 



the requirement is that at each link 

trAB^x) = tiB,j(x) 



(4.5) 



(4.6) 



For the choice A = A (a), the latter equation can be expressed 
in the form 

a cos a + b sina + c = a + c , (4.7) 
where a, b and c are real numbers that depend on B^. Writing 



U,{x) - £/;(*) = A £/„(*) 



S[U„] = S[U^} . (4.3) 



a = \/a 2 + b 2 coscp , b = \/a 2 + b 2 sincp , (4.8) 

the constraint becomes 

cos^cp-a) =coscp. (4.9) 
which is satisfied by the choice 

a = 2cp . 

One also can invent reflections involving the elements 
D% . . .Dm, but this is not necessary for our purposes. 

To illustrate the usefulness of microcanonical reflections, 
we display in Figure 5 the average plaquette value, as a func- 
tion of the number of Monte Carlo sweeps, at several p values 
for both cold and hot starts. In the strong-to-weak coupling 
crossover region, at P = 10, some 4000 thermalization sweeps 
are necessary on a 6 4 lattice volume, when microcanonical 
relections are not employed. Figure 5, right panel, shows the 
improvement if micro-canonical reflections are applied. Ev- 
ery ten Metropolis sweeps through the lattice is followed by 
a series of four reflection sweeps; the thermalization time is 
seen to be greatly reduced. 

Figure 6 is a plot of the plaquette expectation value P(p) = 
(XiPfj V (x)) vs. p, which agrees with the result obtained previ- 
ously in ref. [23]. A rapid crossover is visible near (3=10, but 
according to ref. [23] this is not an actual phase transition. We 
have concentrated our numerical efforts near and just below 
the crossover, in the range P € [9.5, 10]. 

B. Static Quark Potential 

We are interested in demonstrating the "temporary confine- 
ment" property of G(2) gauge theory; i.e. the linearity of the 
static quark potential in some finite distance interval. The 
static quark potential can be calculated from Wilson loops of 
various sizes, but experience tells us that a "thin" Wilson line 




G(2), Metropolis, 6 




P 

FIG. 6: The plaquette expectation values as function of of (3: 6 4 
lattice. 



connecting the quark anti-quark sources has very little over- 
lap with the ground state in quark anti-quark channel. Spa- 
tially smeared links (sometimes called "fat links") are typi- 
cally used for overlap enhancement. APE smearing [24-26] 
adds to particular link a weighted sum of its staples, but the 
subsequent projection onto a SU(N) element usually induces 
non-analyticities. Recently, the so-called "stout link" was in- 
troduced in [27]; in that approach the link variable remains in 
the SU(N) group during smearing. 

All these methods produce time-consuming computer code 



when implemented for the G(2) gauge group. For this rea- 
son, we will use a variant of the smearing procedure outlined 
in [28], in which the spatial links for a given time-slice are 
cooled with respect to a D=3 dimensional lattice gauge ac- 
tion on the time slice. Consider a particular spatial link Ui(x), 
i = 1 ... 3, and apply a trial update U- = G(tt)Ui where 

G(S) = exp{oc a C a } , (4.10) 

and the C a , a = 1 ... 14 are the generators of the G(2) algebra. 
If we denote by M the sum of space-like plaquette variables 
containing the link £/, (x) then the choice 

a a = eti{C a M} , e>0, (4.11) 

will always lead to an increase in trM, for sufficiently small e, 
if the link t/, (x) is replaced by the trial link. Our procedure is 
to start with £«1, and replace links by trial links if trM is in- 
creased by the replacement. Proceeding through the lattice, e 
is tuned so that at least 50% of the changes are accepted. One 
sweep through the lattice is one smearing step. As a stopping 
criterion, one might use the size of AS< 3 ) or a given number 
of smearing steps; our present present results were obtained 
with a fixed number of 70 smearing steps. Let W(r,t) denote 
the expectation value of a timelike, rectangular rxt Wilson 
loop, where r = n s a and t = n t a, and a is the lattice spacing. 
The sides of length r, t are oriented in spacelike and time di- 
rections, respectively. The spacelike links are taken to be the 
smeared links, while the timelike links are unmodified. For a 
fixed value of r = an s , V(r) is determined from a fit of W(r, t), 
computed at various n t , to the form 

c + V(r)a n, « -In (\V(c)\ . (4.12) 

Using smeared links, we find that — In (W(c)) scales linearly 
with n, even for n t as small as n t — 2 (see figure 7 for an ex- 
ample at p = 9.6). The corresponding linear fit yields V(r)a 
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including its error bar. Figure 7 (right panel) shows V(r)a 
as function of rja again for p = 9.6. The line between the 
static quarks (spacelike sides of the Wilson loops) runs par- 
allel to one of the main axis (e.g. (100)) of the cubic lattice. 
In order to check for artifacts arising from the rotation sym- 
metry breaking, the spacelike sides of the Wilson loop were 
also placed along the diagonal axis (110) and (111). From 
Figure 7 it is clear that at this coupling, the rotational symme- 
try breaking effects due to the underlying lattice structure are 
quite small. 

The intermediate string tension a is extracted from a fit of 
the potential V(r) to the form 

V(r) = V Q + or - - , (4.13) 
r 

A sample fit at p = 9.6 is shown in Fig. 7 (right panel). Results 
from 150 independent lattice configurations have been have 
been taken into account. Fitting the data from Wilson loops 
placed along the main axis of the lattice ((100)-data), we find 
(with n = r/a) 

313(5) 

V{r)a = 0.713(7) ■ — + 0.141(2)« , X 2 / dof ~ 2.1 . 

n 

Fitting data from all crystallographic orientations of the Wil- 
son loops, we obtain 

V(r)a = 0.79(1) - ° 34 ^ + 0.13(2)«, x 2 / dof ~9200. 
n 

The large value for % 2 /dof in the latter case arises from the 
facts that (i) only statistically error bars are included in the 
data and that (ii) the error due to rotational symmetry breaking 
is much larger than the statistical one. Both fits are shown in 
figure 7 right panel. 

Our numerical findings provide clear evidence of the linear- 
ity of the potential at intermediate distances r. At large dis- 
tances r, one expects a flattening of the potential due to string 



breaking effects. From experience with SU(2) and SU(3) lat- 
tice theories, we do not expect to observe string breaking from 
measurements of modest-sized Wilson loops. String-breaking 
should be observable in sufficiently large Wilson loops, but 
this would require the rather more sophisticated noise reduc- 
tion methods employed, e.g., in ref. [29]. 

We have calculated the potential with p in the range 
[9.4, 10] on a 12 4 lattice volume. The "raw data" are shown in 
figure 8 (left panel). From a fit of the raw data to eq. (4.13), 
we have extracted the string tension in lattice units, aa 2 , at 
each p value (see table I). In order to avoid any obfuscation 
by rotational symmetry breaking effects, only data from Wil- 
son loops with (lOO)-orientation have been included. The ob- 
tained values for the string tensions can be used to express the 
static potential in physical units, by first subtracting the self- 
energy Vo(P) from the "raw" potential in lattice units, and di- 
viding the result by Vc?a 2 . The result is plotted for several 
P values as function of the distance in units of the physical 
string tension y/ar = n s \J Ga 2 (see figure 8, right panel). A 
fit of the data arising from all P-values to the potential V(r) 
in (4.13) is also included to guide the eye. We observe that 
the data at different P values lies on the same line, and that 
rotational symmetry breaking effects are small. 



C. Finite Temperature Transition 

In SU(2) gauge-Higgs theory there is no unambiguous tran- 
sition between the Higgs phase and the temporary confine- 
ment phase, and no local, gauge-invariant order parameter 
which can distinguish between the two phases (for a recent 
discussion, cf. [30]). In particular, the center symmetry is 
trivial (i.e. Z\), and Polyakov line VEVs are non-zero through- 
out the phase diagram. Nevertheless, there does exist a finite 
temperature transition in this theory, where the Polyakov line 
jumps from a small value to a larger value, and the linear po- 
tential at intermediate distances is lost. 
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G(2), 12 ,(100) (110) (111) 



G(2), 12 ,(100) (110) (111) 





FIG. 8: The static potential in units of the lattice spacing (left panel) and in units of the intermediate) string tension (right panel). 



P 9.5 9.6 9.7 9.8 9.9 10 
"5?" 0.24(2) 0.14(1) 0.102(3) 0.079(1) 0.060(2) 0.047(1) 
a 0.28(5) 0.313(2) 0.318(8) 0.311(2) 0.311(3) 0.309(4) 

TABLE I: Lattice spacing a in units of the string tension o for the values (3; also shown oc(p); 12 4 lattice. 



Denote the Polyakov line as 

P(X) = Y[U4(t,x) 



(4.14) 



The standard order parameter for finite temperature transitions 
is 



P(X) = 77T 



£tr/>(j< 



(4.15) 



on an N 3 x N t lattice, where the temperature is T = 
l/(N t a(fi)). It was observed in refs. [3, 23] that G(2) lat- 
tice gauge theory undergoes a first order transition at a finite 
temperature, indicated by a jump of p(T) at T c . As in SU(2) 
gauge-Higgs theory, this first-order transition cannot be char- 
acterized as a center symmetry-breaking transition. 

Consider a gauge transformation of the Polyakov line 
holonomy 

P'{x) = Q.(x,l)P(x)Q T (x,l) , Q(jf,l) = £l(x,N t ) . 

(4.16) 

By a judicious choice of gauge transformation Q, = TO,, where 
matrices T and £2 are described in Appendix A and Ap- 
pendix C, respectively, we may bring P(x) into the form 



TV P'{x) y f T 1 




U e SU(3) , 



(4.17) 

which emphasizes the SU(3) subgroup. The SU(3) subgroup 
is of interest here because the trace of a Polyakov line is the 
sum of its eigenvalues, the eigenvalues are elements of the 



Cartan subgroup, and the Cartan subgroup of G(2) is the same 
as that of the SU(3) subgroup. We have 



trP(x) = 1 + 2RetrJ7 



(4.18) 



It is therefore sufficient to consider the trace of 3 x 3 subma- 
trix U. If the Polyakov line VEV were exactly zero (the case 
of true confinement, rather than temporary confinement), we 
would have 



1 

- Retrt/ 



1 

6 



(4.19) 



exactly. However, because of color screening, this relation can 
only be approximately true at low temperatures. 
We define 



c(x) = -trU(x) 



(4.20) 



which takes values within the curved triangle shown in fig- 
ure 9. The right panel of figure 9 shows the distribution of 
c(x) for P = 9.7 at low and high temperatures. At p = 9.7 and 
a 6 4 lattice (low temperature case), we clearly observe that 
the points accumulate near the non-trivial center elements of 
SU(3), i.e., 



0.5 



V3 



= -0.5 - i 



V3 



As expected, the scatter plot is symmetric with respect to re- 
flections across the jc-axis. At temperatures T > T c (6 3 x 2 
lattice volume), a strong shift of the data away from the non- 
trivial center elements towards the unit element is clearly ob- 
served. 
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Re c Re c 

FIG. 9: Region of support for values c(3c) in (4.20) for [3 = (left panel). Right panel: Distribution of c(3c) at zero temperature (squares) and 
in the deconfinement phase (circles). 



V. SU(3), Z 3 , OR SU(2) "DOMINANCE"? 

Although the center of the G(2) gauge group is trivial, one 
may ask if, in some gauge, the infrared dynamics is controlled 
by a subgroup of G(2) which does have a center. The motiva- 
tion for this idea goes back to 't Hooft's proposal for abelian 
projection [31]: In an SU(N) gauge theory, imagine impos- 
ing a gauge choice which leaves a remaining U{\) N ~ l gauge 
symmetry. The resulting theory can be thought of as a gauge 
theory with a compact abelian gauge group, whose degrees 
of freedom are "photons" and monopoles, coupled to a set of 
"matter" fields which consist of the gluons which are charged 
with respect to the U(\) N ~ l gauge group. If for some rea- 
son these gluons are very massive, then the idea is that they 
are irrelevant to the infrared physics, which is controlled by 
an effective abelian gauge theory. This is known as "abelian 
dominance". Many numerical investigations have claimed to 
verify the existence of a large mass for the charged gluons. 
Confinement is then supposed to be due to the condensation 
of monopoles in the abelian theory. 

In the same spirit, let us consider fixing the gauge of the 
G(2) theory, in such a way that the action remains invariant 
under an SU(3) subgroup of the G(2) gauge transformations. 
The gauge-fixed theory can be thought of as a theory of SU(3) 
gluons, coupled to vector matter fields (the remaining gluons 
of G(2)) and ghosts. If, for some unknown dynamical reason, 
these vector matter and ghost fields aquire a large effective 
mass, then the linear potential would be largely determined 
by pure SU(3) dynamics, and Z3 vortices associated with the 
center of the SU(3) subgroup would be expected to carry the 
associated magnetic disorder. 

To investigate this idea, we turn to lattice Monte Carlo sim- 
ulations. We have found it useful, for the investigation of 
SU(3) dominance, to use the representation of the G(2) group 
by 7 x 7 complex matrices, as presented in refs. [23, 32]. In 
this representation, any group element Q £ G(2) can be ex- 
pressed as 

g =ZU (5.1) 



where Z is a 7 x 7 unitary matrix which is a function of a 
complex 3-vector K, and 11 is the matrix 

u -{" \) 

where U is a 3 x 3 SU(3) matrix. The construction of the Z 
matrix from a complex 3-vector is described, following ref. 
[23], in Appendix D. The conjecture is that in a suitably cho- 
sen "maximal SU(3) gauge", which brings the Z matrices in 
the g = ZU link decomposition as close as possible, on aver- 
age, to the unit matrix, the information about confinement in 
the intermediate distance regime is carried entirely by the 11 
matrices. This means that the string tension of Wilson loops 
formed, in this gauge, from the 3 x 3 U matrices composing 11 
would approximately reproduce the full intermediate-distance 
string tension, while the string tension of Wilson loops formed 
from the Z variables alone would vanish. If SU(3) dominance 
of this kind is exhibited by G(2) lattice gauge theory, then we 
can make a further reasonable conjecture, which is that the 
confinement information in the U matrices is carried entirely 
by Z3 center vortices. If so, then (i) the Z3 vortices by them- 
selves should also give a good estimate of the full intermediate 
string tension; and (ii) removing vortices from the 11 matrices 
should completely eliminate the intermediate string tension in 
G(2) loop holonomies. 

We define maximal SU(3) gauge to be the gauge which 
minimizes 

* = ll{lli<wi 2 +li:k ; 'wi 2 

x M =l [i=lj=4 i=5j=l 

+ (E + E) l^'WI 2 } (5-3) 

where u'u(x) is the component of the G(2) link vari- 
able (note that Ri = for Z — I7). It is straightforward to 
check that multiplying a link variable on either the left or the 
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right by an SU(3) subgroup element 11 leaves R\ unchanged; 
it is therefore sufficient to use only the g = Z transforma- 
tions in fixing the gauge. This also means that maximizing 
Ri leaves a residual local SU(3) symmetry, so that "maximal 
SU(3) gauge" is a good name for this gauge choice. The next 
step is to use the remaining SU(3) gauge freedom to maximize 



ul\x)+u™{x)+u''{x) 



77/ 



(5.4) 



leaving a remnant local Z3 symmetry in "maximal Z3 gauge." 
This two-step process is reminiscent of fixing to indirect max- 
imal center gauge in SU(N) gauge theory. Details regarding 
our Monte Carlo updating procedure in the complex represen- 
tation of G(2), and the method of fixing to maximal SU(3) and 
maximal Z3 gauges, can be found in Appendix D. 



A. SU(3) and Z 3 Projection 

Having fixed to maximal Z3 gauge, we compute the follow- 
ing observables: 

1. Full Loops. Here we use the full G(2) link vari- 
ables Q = Z 11 to compute Wilson loops, and of course 
gauge-fixing is irrelevant to the loop values. 

2. SU(3) Projected Loops. These are loops computed 
from SU(3) link variables U, extracted from the top 
3x3 block of the 7 x 7 11 link matrices 



U 



11 



(5.5) 



variables Q 1 = Z 1l' . In practice this accomplished sim- 
ply by multiplying the original Q by V\ i.e. construct 
loops from link variables 



(5.7) 



7. SU(3) Vortex-Removed Loops. We construct loops 
from the SU(3)-projected, vortex-removed lattice U 1 = 

z*U. 

Figure 10 shows our results for Creutz ratios of Z3- 
projected configurations at couplings p = 9.5, 9.6, 9.7. The 
straight lines shown are for the intermediate string tension at 
these couplings, extracted from the unprojected lattice and 
listed in Table I above. The agreement of the Z3-projected 
string tension with the full intermediate string tension is quite 
good. 

Z 3 Projected Loops 




FIG. 10: Creutz ratios for Z 3 -projected loops at P = 9.5, 9.6, 9.7. The 
horizontal lines show the values of the corresponding intermediate 
string tensions, computed in section IV on the unprojected lattice. 



U* 



3. Z 3 Projected Loops. These loops are just constructed 
from the z^{x) € Z3 link variables obtained from center- 
projection of the 3 x 3 U -link variables. 

4. Z Loops. These "SU(3)-suppressed" loops are com- 
puted from the Z link variables alone; i.e. we set the 
SU(3) part of the Q = ZU link variables to 11 = 1 7 . 

5. G(2) Vortex-Only Loops. In the decomposition of 
G(2) link variables Q = ZU, replace the SU(3) factor 
11 by the vortex-only element 




(5.6) 



and compute loops with the 7 x 7 link variables Q — 
ZV. 

6. G(2) Vortex-Removed Loops. Here we remove vor- 
tices from the SU(3) elements 11 by the usual procedure 
[33], and use the modified, vortex-removed 11' ', together 
with Z, to construct vortex-removed G(2) lattice link 
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FIG. 11: Creutz ratios for the full, Z3 -projected, G(2) vortex-only, 
and G(2) vortex-removed loops at P = 9.7 

Figure 11 shows the data, at p = 9.7 on a 10 4 lattice, 
for Creutz ratios obtained from full, G(2) vortex-only, G(2) 
vortex-removed, and Z3 projected configurations. Again, de- 
noting the decomposition of the full links in maximal Z3 gauge 
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by g — ZZl, and the Z(3)-projection of u (eq. 5.6) by V, 
then the G(2) vortex-only link variables are g vo = ZV, and 
the G(2) vortex-removed links are g vrem — ZUV^ . It is clear 
from the plot that the G(2) vortex-only data appoaches the data 
obtained from full Wilson loops, while the linear potential is 
absent in the vortex-removed data. 

If the confining disorder resides in the Z3 degrees of free- 
dom, and these in turn lie in the SU(3) projected link vari- 
ables, then SU(3) Dominance is implied. The string tension 
on the unprojected lattice should be seen in the SU(3) projec- 
tion, and removing SU(3) fluctuations from the G(2) link vari- 
ables should send the intermediate string tension to zero. Both 
of these phenomena can be seen in Fig. 12. The suppression 
of the full string tension in the Z loops, which follows from 
suppression of SU(3) fluctuations, is particularly striking. 

Finally, we can remove vortices from the SU(3) link vari- 
ables, and compute Creutz ratios from those lattices. The re- 
sults for p = 9.6 and P = 9.7 are displayed in Fig. 13, and ap- 
pear to be consistent with a vanishing string tension for large 
loops. 

SU(3) Projection at p=9.7 




FIG. 12: Creutz ratios for SU(3) projected lattices, and for G(2) lat- 
tices with the SU(3) link factors set to unity ("SU(3) suppressed"). 
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FIG. 13: Creutz ratios for SU(3)-projected, vortex-removed Wilson 
loops, at p = 9.6, 9.7. 



B. Problems with the Projections 

The success of SU(3) and Z3-projected lattices in reproduc- 
ing the asymptotic string tension of G(2) gauge theory, to- 
gether with the vanishing string tension in "vortex-removed" 
lattices, tends to obscure the fact that these are, after all, 
gauge-dependent observables. That fact doesn't necessarily 
mean that, e.g., the Z3 vortices located via projection are un- 
physical, but neither do the projected results alone prove that 
these are physical objects. The really crucial numerical evi- 
dence in favor of the physical nature of vortices, in SU(2) and 
SU(3) gauge theories, was the correlations that were found 
found between P- vortex location and gauge-invariant observ- 
ables. In particular 

1. Plaquettes at the location of P-vortices have, on the 
unprojected lattice, a plaquette action which is signif- 
icantly higher than average. 

2. A vortex-limited Wilson loop W„(C) is an unprojected 
Wilson loop, evaluated in an ensemble of configurations 
in which n P-vortices pierce the loop of the projected 
lattice. It has been found that as the loop area becomes 
large 



W„(C) 
W (C) 



2nin/N 



(5.8) 



where N is the number of colors. This is the expected 
behavior, if a P-vortex on the projected lattice locates a 
thick vortex on the unprojected lattice. 

3. In an SU(N) gauge theory, a Wilson loop calculated 
in the "vortex-removed" lattice can be equally well de- 
scribed as the expectation value of the product 



W vrem (C) = (Z*(C)Tr[U(C)}) 



(5.9) 



where Z(C), U(C) are loop holonomies on the pro- 
jected and unprojected lattices, respectively. If W vrem 
has a vanishing asymptotic string tension, then this fact 
translates directly into a correlation of the phase of the 
projected loop Z(C), induced by P-vortices, with the 
phase of the gauge-invariant observable Ti[U (C)] on the 
unprojected lattice. 

Unlike the SU(2) and SU(3) cases, we have so far found 
no discernable correlation between gauge-invariant G(2) Wil- 
son loops, and the SU(3) and Z3 projected loops. Plaquettes 
pierced by Z3 vortices have no higher action, on the unpro- 
jected lattice, than other plaquettes, and the VEVs of vortex- 
limited Wilson loops W„ (C) appear to be completely indepen- 
dent of the number of vortices n. Although the string ten- 
sion of vortex-removed loops vanishes in G(2) gauge theory, 
the vortex-removed loops W vrem (C) cannot be expressed as a 
product of vortex loops and gauge-invariant loops, hence the 
VEV of W vrem (C) does not directly correlate the vortex ob- 
servables with gauge-invariant observables. It should be noted 
that in SU(N) gauge theories, vortex removal is a very mini- 
mal disturbance of the lattice, especially at weak couplings. 
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FIG. 14: The G(2) static potential, compared to the static poten- 
tials computed in SU(2)-only and SU(2)-removed lattices, at (3 = 9.6. 
String tensions s = Ga 2 are extracted from fits to eq. (4.13). 



Vortex removal affects the plaquette action only at the loca- 
tion of P-plaquettes, whose density falls exponentially with 
coupling p. In G(2) gauge theory, matters are different. In 
that case Z3 vortex removal is a serious butchery of the origi- 
nal configuration, affecting every plaquette of the lattice. 8 

For these reasons, we think that the results of SU(3) and 
Z3 projection in G(2) gauge theory are misleading. In par- 
ticular, the supposed vortices located by Z3 projection appear 
to not affect any local gauge-invariant observables (i.e. Wil- 
son loops), and on these grounds we believe that the objects 
identified by Z3 projection are simply unphysical. For these 
reasons, we think that the results of SU(3) and Z3 projection 
in G(2) gauge theory are misleading. In particular, the sup- 
posed vortices located by Z3 projection appear to not affect 
any local gauge-invariant observables (i.e. Wilson loops), and 
on these grounds we believe that the objects identified by Z3 
projection are simply unphysical. 9 



In G(2) gauge theory, vortex removal is a matter of multiplying each link by 
a matrix diag[z*/3, 1,2/3] which is not a center element, and which does not 
commute with the factor we have denoted by Z . In general this procedure 
will change the plaquette action, in an uncontrolled way, at every plaquette 
on the lattice. This is in complete contrast to SU(N) gauge theories, where 
multiplication of link variables by projected center elements looks locally, 
almost everywhere, like a Zv gauge transformation, and changes plaquette 
actions only at the location of P-plaquettes. 
9 A possibly related issue has been studied recently by Pepe and Wiese [34]. 
These authors simulate G(2) lattice gauge-Higgs theory, where the Higgs 
field can spontaneously break the symmetry down to SU(3). An expec- 
tation value for the Higgs field gives a mass to six of the fourteen G(2) 
gluons, while eight other gluons, associated with the SU(3) subgroup, re- 
main massless. The masses of the massive gluons can be adjusted, from 
zero to infinity, by the hopping parameter in the Higgs Lagrangian, and 



One might ask if the reason for the ultimate failure of the 
SU(3) and Z3 projections lies in the choice of the subgroup. 
If we would calculate Wilson loops from the 7x7 SU(3)- 
projected link variables U, rather than the 3x3 SU(3) link 
variables U, then those Wilson loops would not have an area 
law, due to the unit element in the 7x7 W matrix; i.e. the 
existence of an element of the fundamental representation of 
G(2) which transforms as a singlet under the SU(3) subgroup. 
Suppose that we instead look for the smallest subgroup of 
G(2) such that all elements of the fundamental representa- 
tion of G(2) transform as non-singlets under that subgroup. 
It turns out that the SU(2) subgroup defined by the generators 
C$,Cg,C\o (see Appendix A) satisfies this condition; the fun- 
damental representation of G(2) transforms, under this sub- 
group, as a triplet and two doublets. 

We then impose a "maximal SU(2)" gauge, and calculate 
the Wilson loops of the SU(2)-only and SU(2)-removed link 
variables. The results are impressive, and are displayed in Fig. 
14. However, up to now we have not found any significant 
correlation between SU(2)-only and G(2) Wilson loops, and 
in this respect the SU(2) projection appears to be also prob- 
lematic. 

We believe that the reason that SU(3), Z 3 , and SU(2)- 
projected loops apparently fail to correlate with gauge- 
invariant observables is that these projections have no gauge- 
invariant motivation from the beginning. We may reasonably 
suppose that a gauge theory with gauge group G is "domi- 
nated" by the degrees of freedom associated with a subgroup 
H, if, independent of any gauge-fixing, the gauge-invariant 
string tension depends only on the representation of the sub- 
group. This is the case for the Zn center subgroup of SU(N) 
gauge theories, at large distance scales. Center projection is 
then a method for isolating the large-scale confining fluctu- 
ations. In contrast, the SU(3), Z3, and SU(2) subgroups of 
G(2) have no relation to G(2) Casimir scaling, at intermediate 
distances, nor are any of these subgroups singled out at large 
scales, where the string tension vanishes for any representa- 
tion. 10 In fact, similar remarks can be made regarding abelian 
projection in pure SU(N) gauge theories, which we used to 
motivate the subgroup projections of this section. The abelian 
U(1) A '~ 1 subgroup cannot easily account for Casimir scaling 
at intermediate distances, while at asymptotic distances the 
string tension of an abelian-projected loop in fact depends 
only on the N-ality, rather than the abelian charge of the loop 
[35]. 



in this way the gauge-Higgs theory can interpolate between pure G(2) and 
pure SU(3) gauge theory. The question addressed in ref. [34] is whether 
the deconfinement transition of SU(3) gauge theory, which is certainly a 
Z3 symmetry-breaking transition, is continuously connected, in the phase 
diagram, to the first-order finite-temperature transition in G(2) gauge the- 
ory. The answer appears to be no, and the authors argue that in fact the two 
transitions have different origins. This conclusion seems consistent with 
our finding that observables defined in SU(3) and Z3 projections appear to 
be physically irrelevant in pure G(2) gauge theory. 

The situation would be different in a G(2) gauge-Higgs theory, with G(2) 
broken to SU(3). 
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VI. CONCLUSIONS 

In this article we have presented a unified picture of vac- 
uum fluctuations in G(2) and SU(N) gauge theories, in which 
a two-dimensional surface slice of the four-dimensional lattice 
can be decomposed into domains corresponding to elements 
of the gauge group center. Domains associated with center 
elements different from the identity are just two-dimensional 
cross-sections of the usual Zn vortices; the new feature is that 
we also allow for domains associated with the unit center el- 
ement, which exists in G(2) as well as SU(N). This domain 
structure accounts for the representation dependence of the 
asymptotic string tension on the center subgroup, which in 
the case of G(2) implies a vanishing asymptotic string ten- 
sion. The existence and Casimir scaling of the linear potential 
at intermediate distances is explained by random spatial vari- 
ations of the color magnetic flux in the interior of the center 
domains. 

The picture we have outlined is a kind of phenomenology. 
Domain structure is arrived at by simply asking what features 
of typical vacuum fluctuations, in a Euclidean 4-volume, can 
account for both Casimir scaling of intermediate string ten- 
sions, and color-screening in the asymptotic string tensions. 
Of course one would like to support this picture by either 
analytical methods or numerical simulations. There is, in 
fact, abundant numerical evidence for the relevance of Z2 and 
Z3 center vortices in SU(2) and SU(3) lattice gauge theories 
[1, 2], but unfortunately the usual tools which are used to lo- 
cate such vortices in lattice simulations (i.e. center projection 
in maximal center gauge) are not of much use in locating do- 
mains associated with the unit element. At the moment we can 
only appeal to a number of analytical studies of the one-loop 
effective action of vortex configurations, which show that this 
action is stationary for magnetic flux quantized in units cor- 
responding to the gauge group center elements, including the 
unit element [20-22]. Thus the dynamics underlying center 
vortex formation may ultimately be the same for the unit and 
non-unit element varieties. The picture we have presented for 
Casimir scaling at intermediate distance scales is motivated by 
the notion of dimensional reduction, and by various approx- 
imate treatments of the Yang-Mills vacuum wavefunctional 
[7,11,13,15-17]. 

We have also calculated certain static properties of G(2) lat- 
tice gauge theory, namely, the intermediate string tension, and 
the distribution of Polyakov line holonomies at zero and finite 
temperature. We see very clear evidence of a linear potential 
in the fundamental representation, but we have not attempted 
to compute the potential for higher representations, needed to 
verify Casimir scaling, or the operator-mixings needed to see 
string breaking. These are left for future investigations. The 
eigenvalues of the Polyakov line are those of the SU(3) sub- 
group, and we have shown that the finite temperature phase 
transition is accompanied by a shift in the concentration of 
the Polyakov line eigenvalues from non-trivial to trivial SU(3) 
center elements. 

Finally, we have investigated projection of G(2) link vari- 



ables to various subgroups. These projections, to SU(3), Z3, 
and SU(2) subgroups, so far lack the gauge-invariant motiva- 
tion that exists for center projection in SU(N) gauge theories, 
where the representation dependence of the asymptotic string 
tension is given entirely by the representation of the Zn sub- 
group. The SU(3), Z3, and SU(2) subgroups of G(2) have 
no natural connection to either G(2) Casimir scaling at inter- 
mediate distances, where a non-zero string tension exists, or 
to the asymptotic regime, where the string tension vanishes. 
The fact that we have found no correlation whatever, between 
the projected Wilson loops and gauge-invariant observables 
in G(2) gauge theory, probably reflects the fact that these sub- 
groups do not appear naturally in any gauge-invariant set of 
observables. The only subgroup which really does appear to 
be relevant to pure non-abelian gauge theories is the center 
subgroup, and this relevance persists even in cases such as 
G(2) gauge theory, where the center subgroup consists of only 
a single unit element. 
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APPENDIX A: G(2) ALGEBRA 

In the following, we will study the fundamental representa- 
tion of G(2) in more detail. In this representation, the genera- 
tors can be chosen as real 7x7 matrices. The group elements 
U satisfy the constraints 

UU T = 1 , T abc = T def U da U eb U fc , (Al) 

where T a t c is a total anti-symmetric tensor [3], The group 
elements are generated by the Lie algebra: 

U = exp j-i £ p a /C fl | , p e %. u . (A2) 

An explicit realization of the generators iC a was presented 
in [36]. For the readers convenience, these matrices are pro- 
vided below. For this choice of representation, these elements 
are 

Tm = T\45 = Tne = T 2 46 = T251 = T^i = r 365 = 1 . 

(A3) 

All linear combinations of the generators C a form the vector 
space L of the Lie algebra. 
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The embedding of G(2) in the group SO(7) is generated by the 14 Lie algebra elements Ck,k— 1 , . . . , 14: 



Ci = 



C 3 = 



C 5 



C 9 



Ci 3 = 



/o 

















o \ 






























































-1 

















-1 

















1 








\o 








1 








o / 


/° 

















\ 
























































-1 

















1 





























-1 


\o 














1 


o / 



c 7 = 



/0 









1 

\0 

/0 





1 

1 



\0 

/° 

2 










-1 



















\ 

-1 



1 

7! 





-1 







-2 



Cn = -t= 



1 

V3 



1 

V3 








Vo 

fO 





2 





Vo i 

/0 





1 

1 

2 

\0 










-1 













-1 

-2 









1 





0/ 

0\ 











0/ 









-1 






c 2 



c 4 







1 











°\ 




1 




0/ 

\ 

-1 












-1 










-1 









-2 









/ 
°\ 







0/ 



c 6 = 



c 8 = 



/0 











\0 

/° 







Vo 

/° 







Vo 



1 

71 









-1 











1 





1 








-1 















0\ 















1 















1 









































-1 









-1 








J 















0\ 















1 












-1 












1 












-1 








































J 






f° 














\ 








-2 











2 






















1 
















-10 



















-1 




Vo o 





1 





/ 



Cio 



1 

7! 



Cl2 = vf 



Cu — 



V3 



/0 


-2 


\ 














2 

















-1 














-1 








10 





Vo 





10 


/ 


/° 





0-2 


0\ 











1 











1 














2 














-1 








Vo 





-10 


0/ 


f° 








- 2 \ 








10 











0-10 








-1 














10 

















V2 








/ 



(A4) 



The two matrices C3 and Cs generate the Cartan subgroup of G(2). Furthermore, there are 6 St/ (2) subgroups generated by the 
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elements: 



C\, C 2 , Y\ ■= C3, 

C 4 , C 5 , y 2 :=7(\/3C 8 + C 3 ), 

C 6 , C 7 , y 3 :-3(-\/3C8 + C 3 ), 

V3C9, \/3Cio, I4 := V3C8, 

v^Cn, 73Ci2, y 5 :-3(-V3C 8 + 3C 3 ), 

V3Ci 3 , V3Ci4, y 6 :=Hv / 3C 8 + 3C 3 ), 



[Ci,C2]-2Fi, 
[C 4 ,C 5 ] -2F 2 
[C 6 ,C 7 ] =2T 3 
[V3C 9 ,\/3Cio] -2F 4 
[^Cii,^Ci2]=27 5 
[V3C 13 ,V3C H }=2Y 6 



(A5) 



The first three SU(2) subgroups form 4 dimensional real representations of St/ (2) as subgroups of 5(9(7) and they generate an 
SU(3) subgroup of G(2). The representations of the remaining three SU(2) subgroups (as subgroups of 5(9(7)) are a little bit 
more complicated - they are 7 dimensional but they are reducible: the representation space is the direct sum of a 3 dimensional 
adjoint representation and a 4 dimensional real representation (as in the case of the first three SU(2) subgroups). 
The various SU (2) subgroups of G(2) C 5(9(7) are obtained by exponentiation: 
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The elements of the two remaining SU (2) subgroups look quite similar. 

Alternatively, one may equally well use a complex representation of G(2), as in refs. [3, 23] and in section V of this article. 
The two representations can be related by a similarity transformation Q = VC^V^, where V is the unitary matrix 
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and in particular 



Ci - 



c 4 



c 7 




,c 3 = 



1 Q 





0-10 

10 




















— 2/03 




























/03 




(A10) 



Finally one can use the (unitary) permutation matrix 



/I 000000X 

1 

1 

1 

1 

1 

\o 1 0/ 



(All) 



to rearrange the rows and columns of C\,...,C% such that it is obvious that the decomposition of the 7-dimensional St/ (3) 
representation is given by {1} {3} ® {3}. 
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APPENDIX B: EULER DECOMPOSITION OF G(2) 



In order to implement a Metropolis update, we will use prototypes of G(2) group elements to maneuver through group space. 
These are: 
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which are "pure" SU(2) subgroup elements. Furthermore we use: 
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Following [36], the Euler angle representation of an arbitrary G(2) element G is given by: 



G = Da(ai)D g (a2)D & (a 3 )D 3 (a4)D2(a 5 )D 3 (a 6 ) 

Dn(a 7 )£)5(oc8) £>3(a9)£> 2 (aio)£>3(aii) £>8(ai2)£> 9 (ai 3 )Z)8(ai4) > 
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where the range of parameters is found to be: 

a\ , A4, «9, an G [0,2jc[, «2, «5, aio, «13 € [O^ti a 3, «6, an, «14 G [0,jc[ . (B4) 

and 

a% G [0,J[, a 7 G [0,|flg[. (B5) 



APPENDIX C: POLYAKOV LINE AND THE CARTAN 
SUBGROUP 

1. Higgs field and unitary gauge 

As pointed out by Holland et al. [3], a fundamental Higgs 
field will break G(2) gauge symmetry to SU(3). This can be 
easily anticipated resorting to the real representation of G(2). 
Also for later use, we will here provide the sequence of 5(9(2) 
rotations Dj(a) with which the 7 component Higgs field can 
be rotated into the 1 -direction: 



where 



* 

v*y 



Dio(ai) 
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(CI) 



Here, the symbol * denotes an arbitrary real number. The dif- 
ferent angles a, can be easily computed step by step. To be 
precise, we here provide some details. Let us consider the vec- 
tor (x,y) T the y component of which we would like to rotate 
to zero. There are two possibilities for that: After the rotation 
of the vector, its x component can be either positive or nega- 
tive. Throughout this section, we adopt the "minimal" choice, 
which preserves the sign of the x component: 



cos(j) — sin(j) 
sincj) coscj) I \y 



sign(x) \Jx 2 +y 2 



COS(j) = 



\fx 2 +y 2 ' 



sin(j) = 



-sign(x) y 



Now we will show that a non-zero Higgs field breaks the 
symmetry down to SU(3). To this end we have to calculate the 
stabilizer of the vector (j)o = (l,0,0,0,0,0,0) r , i.e. we have to 
determine the G(2) group elements g fulfilling 



t»o = <>o - 



(C2) 



If we are only interested in continuous subgroups as stabilizer 
we only need to analyze the neighborhood of the identity. In- 
finitesimally eq. (C2) translates to 

= £p fc Q(h = ^=(0,p9,Pio,Pn,Pi2,Pi3,Pi4) r , (C3) 
k V3 

i.e. Pg = pio = Pn = P12 = Pi3 = Pi4 = 0. But this means that 
the only continuous subgroup of G(2) leaving (j)o invariant is 
generated by C\ ,k — 1, ... ,8, i.e. it is the SU(3) subgroup. 
With our considerations we cannot exclude other elements of 
G(2) (forming a discrete subgroup) which fulfil eq. (C2). 



2. The Cartan subgroup 

The Polyakov line P homogeneously transforms under the 
gauge transformation £1. Let H be an element of G(2) with 



Q.HQ. 1 



(C4) 



For an arbitrary element P of G(2), we introduce a correspond- 
ing 7-dimensional vector with the help of the constraint con- 
stants Tabc (A3) by 

(P)a := Tabc Pbc ■ (C5) 

Rewriting (Al) as 

Tauv Qub &vc — Qam Tmbc i 

we easily show that the vectors p, h, constructed from P and 
H respectively, are related by 



Slh. 



(C6) 



In the last subsection, we have already verified that any (real) 
7-dimensional vector can be rotated to e\ = (l,0,0,0,0,0,0) r 
direction using Q. G G(2) only. 
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After a suitable choice of gauge, we now consider Polyakov 
lines P £ G(2) which satisfy p°^e\. Using the constraint 
(Al), we find: 

Tdef Pda Peb T^abc Pfc > 

and in particular (after summing over b = f) 
P T p = p. 

Hence, all elements G £ G(2) with g°^e\ constitute the sub- 
group of rotations which leave the vector e\ invariant. For 
£^3 C SU(3), where St/ (3) is the subgroup of G(2) spanned 
by the generators C\ . . .C%, we find 

e\ — Q.3 e\ , 

implying that these rotations are identified with the SU(3) C 
G(2) subgroup. 

Finally, there are elements £^3 £ SU (3) C G(2) which trans- 
form P £ SU (3) to its center elements. These elements con- 
stitute the rank 2 Cartan subgroup of G(2). 

APPENDIX D: COMPLEX REPRESENTATION AND 
GAUGE-FIXING 

Matrix Z has the form 

/ C pK D* \ 

z = l j£ -pK T (Di) 

\ D pK* C* J 

where x =|| K || 2 , p = V2/(l +x), and C,D are the 3 x 3 ma- 
trices 

c = l{ h -wU} 

W S 

° = -A~A2 (D2) 

with 

M = KK^ , S = K*K f 

W = £ a fryKy , A= Vl +x (D3) 

Each group element is specified by fourteen parameters: eight 
for the SU(3) matrix U, and another six for the complex 3- 
vector K. 

Our Monte Carlo updates combine a Cabbibo-Marinari up- 
date of the link variables, using group elements in the SU(3) 
subgroup of G(2) (i.e. the 11 matrices), followed by a gauge 
transformation by the Z matrices, chosen at random from a 
lookup table. 1 1 The lookup table of several thousand Z matri- 
ces is generated stochastically, by choosing the real and imag- 
inary parts of each component of the complex K vector from a 



1 1 The idea of combining SU(3) Cabbibo-Marinari updates with random G(2) 
gauge transformations was suggested to us by M. Pepe. 



uniform distribution of random numbers in the range [—1,1]. 
For each Z matrix entered into the table, one also enters its in- 
verse, Z~ l = Z 1 . To carry out the Cabbibo-Marinari updates, 
it is necessary to generate a matrix 




(D4) 



where the 3 x 3 A matrix belongs to one of the three standard 
SU(2) subgroups of SU(3), and has non-zero elements only in 
the i, j'-th rows and columns . Let Qi be the link variable at 
link Z, F the associated sum of staples, and R = Q[F. Also 
denote ;' = ;' + 4, / = j + 4, and introduce the 2 x 2 matrices 

V a j' a jj J 

r = ( R "+ R h' R >> + R h> ) 

V R .i i + R fi> R jj +R *ff ) 

3 

17=1 

3 

r = rol2-i£r„a„ (D5) 

«=i 

where the a„'s are real, and the r n 's are, in general, complex. 
Then 

B 2B 3 

^ReTr[« J R] = -^^a„Re(r„) (D6) 

' ' n=0 

From here, one generates a n by the usual heat bath. This is 
done, at each site, for each of the three usual SU(2) subgroups 
ofSU(3). 

We have checked that the plot of plaquette energies vs. cou- 
pling B, generated by this procedure, agrees with the corre- 
sponding result arrived at by a more conventional Metropolis 
updating method, described in section IV above. The plot also 
agrees with the curve displayed in Fig. 2 of ref. [23]. 

To fix to maximal SU(3) gauge, a simple steepest-descent 
algorithm seems to be adequate. We begin with a number of 
"simulated annealing" sweeps at zero temperature, i.e. at each 
site a random complex 3 -vector K is generated, and the cor- 
responding Z is used as a trial gauge transformation. If R\ is 
lowered, the change is accepted. This is followed by steepest 
descent sweeps. Denoting the components of the K vector as 
Kj = Xi + k I+ 3, we numerically compute the gradient, at any 
given site 

dR 

(VR) i = jr i = 1-6 (D7) 
axi 

Then we set 

Xi = -(VR)ie (D8) 

where e is gradually reduced, as the iterations proceed, so as 
not to overshoot the minimum too often. Although conjugate 
gradient methods are preferable to steepest descent, in this 
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case the number of iterations required to converge to a local 
minimum of R± are not excessive. The main cost in gauge fix- 
ing comes at the second step, i.e. fixing to maximal Z3 gauge. 
We have found that steepest descent is inadequate for max- 
imizing R2- Instead, at each site, we use a standard (quasi- 
Newton) optimization routine to maximize R2 with respect to 
the eight Euler angles [37] specifying an SU(3) gauge trans- 
formation at that site, and proceed to fix the gauge by ordinary 
relaxation. 

In order to do the SU(3) projection after maximal Z3 gauge 
fixing, it is necessary to factor each link variable, which is a 
7x7 G(2) matrix, into the product Z 11 . This is carried out as 
follows: For a given G(2) matrix g , define the matrix 

M(K) =z\K)g (D9) 

Then there is some choice of K such that M(K) = 11, and 
M(K) has the block-diagonal form shown in eq. (5.2). In that 
case there are 30 matrix elements of M(K) which vanish, and 
one that equals unity, but there are only six independent de- 
grees of freedom in K. It is sufficient to compute K from only 
six of the block-diagonal conditions. We have chosen these, 



somewhat arbitrarily, to be 



MK) 


= Re[M 44 


-1]=0 


fi(K) 


= Re[M 35 ] 


= 


MK) 


= Re[M 26 ] 


= 


MK) 


= Re[M 52 ] 


= 


MK) 


= Re[M 61 ] 


= 


MK) 


= Re[M 73 ] 


= 



Solving the set of equations /, ■ = determines K. This set can 
be solved numerically; it can be checked that the correspond- 
ing 11 = M(K) satisfies all of the remaining block-diagonal 
conditions, and the 3 x 3 U matrix and its complex conjugate 
in the non-zero blocks are unitary. 

The final step, projection of the 3 x 3 SU(3) matrix U to 
its nearest center element, is standard. The prescription is to 
select the SU(3) center element z G Z3 such that zJ.3 is closest 
to the 3 x 3 SU(3) matrix U ; i.e. choose the z which maximizes 
Re[z*Tr[f/]]. 



[1] J. Greensite, Prog. Part. Nucl. Phys. 51, 1 (2003) [arXiv:hep- 
lat/0301023]. 

[2] M. Engelhardt, Nucl. Phys. Proc. Suppl. 140, 92 (2005) 

[arXiv:hep-lat/0409023]. 
[3] K. Holland, P. Minkowski, M. Pepe, and U.-J. Wiese, Nucl. 

Phys. Proc. Suppl. 119, 652 (2003) [arXiv:hep-lat/0209093]; 

Nucl. Phys. B668, 207 (2003) [arXiv:hep-lat/0302023]. 
[4] K. Fredenhagen and M. Marcu, Phys. Rev. Lett. 56, 223 (1986). 
[5] L. Del Debbio, M. Faber, J. Greensite and S. Olejmk, Phys. 

Rev. D 53, 5891 (1996) [arXiv:hep-lat/9510028]. 
[6] M. Faber, J. Greensite and S. Olejmk, Phys. Rev. D 57, 2603 

(1998) [arXiv:hep-lat/9710039], 
[7] J. Greensite, Nucl. Phys. B 158, 469 (1979). 
[8] P. Olesen, Nucl. Phys. B 200, 381 (1982). 
[9] J. Ambj0rn, P. Olesen, and C. Peterson, Nucl. Phys. B240, 189 

(1984). 

[10] M. B. Halpern, Phys. Rev. D 19, 517 (1979). 
[11] P. Mansfield, Nucl. Phys. B 418, 113 (1994) [arXiv:hep- 
th/9308116]. 

[12] A. I. Shoshi, F. D. Steffen, H. G. Dosch and H. J. Pirner, Phys. 

Rev. D 68, 074004 (2003) [arXiv:hep-ph/02 11287]. 
[13] J. Greensite, Nucl. Phys. B166, 113 (1980). 
[14] J. Greensite and J. Iwasaki, Phys. Lett. B 223, 207 (1989). 
[15] D. Karabali, C. J. Kim and V. P. Nair, Phys. Lett. B 434, 103 

(1998) [arXiv:hep-th/9804132]. 
[16] R. G. Leigh, D. Minic and A. Yelnikov, Phys. Rev. Lett. 96, 

222001 (2006) [arXiv:hep-th/0512111]. 
[17] L. Freidel, arXiv:hep-th/0604185. 

[18] M. Faber, J. Greensite and S. Olejmk, JHEP 0006, 041 (2000) 

[arXiv:hep-lat/0005017]. 
[19] J. Greensite and S. Olejmk, JHEP 0209, 039 (2002) [arXiv:hep- 

lat/0209088]. 

[20] J. D. Lange, M. Engelhardt and H. Reinhardt, Phys. Rev. D 68, 



025001 (2003) [arXiv:hep-th/0301252]. 
[21] D. Diakonov and M. Maul, Phys. Rev. D 66, 096004 (2002) 

[arXiv:hep-lat/0204012] 
[22] M. Bordag, Phys. Rev. D 67, 065001 (2003) [arXiv:hep- 

th/0211080]. 

[23] M. Pepe, PoS LAT2005, 017 (2006) [Nucl. Phys. Proc. Suppl. 

153, 207 (2006)] [arXiv:hep-lat/05 10013], 
[24] M. Albanese et al. [APE Collaboration], Phys. Lett. B 192, 163 

(1987). 

[25] M. Teper, Phys. Lett. B 183, 345 (1987). 
[26] S. Perantonis, A. Huntley and C. Michael, Nucl. Phys. B 326, 
544 (1989). 

[27] C. Morningstar and M. J. Peardon, Phys. Rev. D 69, 054501 

(2004) [arXiv:hep-lat/0311018]. 
[28] K. Langfeld, Phys. Rev. D 69, 014503 (2004) [arXiv:hep- 

lat/0307030]. 

[29] S. Kratochvila and P. de Forcrand, Nucl. Phys. B 671, 103 

(2003) [arXiv:hep-lat/0306011]. 
[30] J. Greensite and S. Olejmk, Phys. Rev. D 74, 014502 (2006) 

[arXiv:hep-lat/0603024] . 
[31] G. 't Hooft, Nucl. Phys. B153 (1979) 141. 
[32] A. J. Macfarlane, Int. J. Mod. Phys. A17 (2002) 2595. 
[33] P. de Forcrand and M. D'Elia, Phys. Rev. Lett. 82 (1999) 4582 

[arXiv: hep-lat/9901020]. 
[34] M. Pepe and U.-J. Wiese, arXiv: hep-lat/06 10076. 
[35] J. Ambj0rn, J. Giedt and J. Greensite, JHEP 0002, 033 (2000) 

[arXiv:hep-lat/9907021]. 
[36] S. L. Cacciatori, B. L. Cerchiai, A. Delia Vedova, G. Ortenzi 

and A. Scotti, J. Math. Phys. 46, 083512 (2005) [arXiv:hep- 

th/0503106]. 
[37] M. Byrd, arXiv:physics/9708015. 



